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Wo first determine the AiV 5-wave phase shifts so as to reproduce the experimental A separation 
energies of A = 3, 4 A-hypernuclei (^H, A H, ^He, ^H*, and A He*), and then construct three phase- 
equivalent AN potentials with different central repulsion. By the stochastic variational method with 
correlated Gaussian basis we perform an extensive calculation of ab initio type for the hypernuclei of 
' up to A— 6. The binding energies and the sizes of the A-hypernuclei are very insensitive to the type 

of the phase-equivalent AiV potentials. We use two different AA potentials which both reproduce 
A_Baa(aaHb) reasonably well. Any combination of these AiV and AA potentials predicts hitherto 
undiscovered particle-stable bound states, aaH, aaH and A^He: Predicted values of -Baa are about 
^ ■ 0.4, 5.5 and 6.3 MeV, respectively. The binding energy of aaH is so small that its possibility crucially 

Q depends on the strength of the AA interaction. The binding energies of both A He and A^He are 
calculated to be strongly overbound compared to experiment. In relation to this well-known anomaly 
£ — ■ we examine the effect of the quark substructure of N and A on their binding energies. The effect is 

£\J ' negligible if the baryon size in which three quarks are confined is smaller than 0.6 fm, but becomes 

appreciable, particularly in AAHe, if the size is taken to be as large as 0.7 fm. We discuss the extent 
to which the nucleon subsystem in the hypernuclei changes by the addition of A particles. The 
charge symmetry breaking of the AN potential is phenomenologically determined and concluded to 
be weakly spin-dependent. 
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I. INTRODUCTION 



Strangeness or hypernuclear physics has been attracting increasing interest in recent years Q . Experimental progress 
on both a production of hypernuclei with strangeness S = — 1 and —2 and an analysis of their decay modes will bring 
rich information on the baryon dynamics involving strange particles. It may be particularly interesting to make clear 
to what extent the properties of a strange particle in nuclear medium change from those in a free space. 
^ , The knowledge of the hyperon-nucleon (YN) interaction is crucially important for a further development of the 
hypernuclear physics. Because of the limited data on A-nucleon (AiV) scatterings, no reliable phase shift analysis 
has been performed yet. Even the cr^-cr^ part of the AiV interaction has not been determined precisely. Though 
one-boson-exchange (OBE) models [plpl or quark cluster models |I| for YN interactions have produced useful results 
on the interactions and have shown some predictive power, the behavior of the 1 Sq or the 3 <Si phase shift is not always 
the same in these models. At present it is not possible to pin down realistic potentials which can be used in ab initio 
hypernuclear structure calculations. 

In view of this circumstance, hypernuclear spectroscopic study will be aided by modeling effective AiV and AA 
potentials. One may determine the potentials phenomenologically so as to reproduce the binding energies of light 
A- and AA-hypernuclei before theoretically sound, realistic potentials become available. In fact, the early work || 
performed more than a quarter of a century ago with the effective AiV potential demonstrated that the observed 
binding energy of ^He is anomalously small. This problem still remains an enigma Q despite of several attempts at 
resolving this anomaly by the mechanism of e.g., tensor forces or AN-HN couplings. Since the excited states of the 
A = 4 isodoublets (^H* and y\He*) had not been well established in those days, it will be worthwhile to redetermine 
effective AA?" potentials which include the charge symmetry breaking (CSB) effect. 

The information on effective AiV interactions becomes available together with the progress in technique for a precise 
solution of few-body systems. A detailed, systematic analysis of the binding mechanism of light hypernuclei may make 
it possible to determine the S'-wave strength of the AiV interaction. As is well-known, the binding energies of the 
A = 3, 4 hypernuclei strongly depend on the relative strength of attraction between 1 So and 3 Si states. If the 1 Sq 
attraction is only slightly stronger than the 3 5i one as in the case of the Nijmegen model F ||(NF) or Dalitz et al. 

H potential, we would have not only the ground state ( J T = \ ) but also a particle-stable excited state ( J n — § ) 
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of ^H. The latter has not experimentally been confirmed yet. This type of potentials predicts a too small energy 
difference of the two states (0 + and 1 + ) for the A = 4 isodoublets. On the contrary, one of the versions of the quark 
cluster model, FSS j|, gives a very strong attraction in 1 5o and leads to the overbinding of A H as well as a too large 
energy difference for the two states in the A = 4 system. Recently the simplest three baryon system YNN with 
S = — 1 has been studied Q in the Faddeev method, which shows that only the Nijmegen soft core (NSC89) potential 
among others, reproduces the experimental J3a(^H) value. The work has also shown clearly the importance for 
the AiV-EiV conversion and its 3 5i- 3 £>i tensor coupling for producing the bound state of |H. The new version of the 
Nijmegen soft core model (NSC97 has been constructed on the basis of these theoretical developments. 

The purpose of the present study is threefold: First we determine phenomenologically such effective A7V central 
potentials that reproduce the binding energies of s-shell hypernuclei (%H, ^H, ^He, \H*, ^He*). Secondly, we explore 
the possibility of whether AA-hypernuclei with A = 4, 5 form bound states or not by employing the available AA 
potentials. This will be useful for experimental search for these AA-hypernuclei. A confirmation of the lightest AA- 
hypernuclei as well as the H-particle will give us useful information on the strength of the AA interaction. Thirdly, 
we investigate to what extent the quark substructure of the baryons plays a role in reducing the A separation energy 
of the hypernuclei with A > 5. 

The third point mentioned above is expected to reveal, through hypernuclear studies, a unique role of the quark 
substructure of baryons. The binding energies of the A-hypernuclei have so far been calculated on the premise that 
they ought to be understood provided two- and three-body baryonic interactions are precisely known. This seems 
quite natural because the A-particle can be distinguished from the nucleon at the baryon level. However, if we take 
into account the quark substructure of the baryon, the Pauli principle acting at the underlying quark level produces 
certain constraints on the dynamical motion of the hypernuclear system, because both N and A contain u and d 
quarks. This quark Pauli effect has recently been studied in Refs. jll[| and |Q as a function of the baryon size b in 
which three quarks are confined. The reduction of the binding energy of ^He and ^^He is found to occur at b > 0.6 
fm. 

It is of vital importance for achieving the above objectives that we can obtain a precise solution of few-body 
problems. A fair description of the correlated motion is an essential ingredient to get accurate solutions. As tested 
in a variety of examples in atomic and nuclear few-body problems [OUlJ], the correlated Gaussian (CG) basis leads 
to a virtually exact solution. This basis explicitly describes the correlated motion of the constituent particles, and its 
simple form makes it possible to obtain matrix elements analytically. We therefore use a variational trial function given 
as a combination of the CGs, and increase the basis dimension one by one till a practical convergence is attained. The 
optimization of the nonlinear parameters contained in the basis functions is performed by selecting the best among a 
number of candidates randomly chosen. This gambling procedure is called the stochastic variational method (SVM) 

us 

We will focus on central components, both spin- independent and spin-dependent, of the potential, so the systems 
considered here are limited to those hypernuclei whose total orbital angular momentum is well approximated by L = 0. 
We presume that the effects of tensor forces, three-body forces and A-S conversion etc. are effectively replaced by 
a central force in fitting the binding energies of the s-shell hypernuclei. The hypernuclei with A < 6 are treated, 
without any compromise, as an A-baryon system and their energies are calculated by fully taking into account the 
dynamics of the constituent baryons. 

We present in sect. 2 the formalism of the present study. The effective A7V potential is determined in subsect. 2.1 by 
starting from potentials phase equivalent to the Nijmegen model || or the quark cluster model Q. The AA potentials 
used in the present study are briefly mentioned in this subsection. The variational trial function with the CGs is given 
in subsect. 2.2. The optimization of the parameters of the trial function is carried out as described in subsect. 2.3. 
Results of calculation are given in sect. 3. The A-hypernuclei with A < 5 and AA-hypernuclei with A < 6 are discussed 
in subsects. 3.1 and 3.2, respectively. The A density distribution and the AA-correlation function are calculated in 
subsect. 3.3 in order to exhibit the structure of the hypernuclei. We also analyze the extent to which the nuclear core 
changes by adding A particles. The CSB term of the AN potential is discussed in subsect. 3.4 and its strength is 
determined phenomenologically. We consider in subsect. 4.1 a special configuration constrained by the Pauli principle 
at the quark level. Results for ^He, and ^He are presented in subsect. 4.2. Summary is given in sect. 5. Some 
formulas for the matrix elements of the CGs are collected in appendix. 

II. THE FORMALISM 
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A. Potentials 



As mentioned in the introduction, our aim in this subsection is to determine effective AJV central forces that 
presumably include various effects such as the tensor force, the three-body force and so on. The Minnesota potential 
fig] is such a type of central potential widely employed as the NN potential. The Minnesota potential is given in a 
form 

V = jv* + \{l + P°)Vt + \{l - P°)V S } \\u+\{2- u)P r } , (1) 
V R = V 0R e~^ r2 , V t = ~V ot e- K ^ 2 , V s = -V 0s e~^ , 

where P a and P r are spin- and space-exchange operators, respectively. An adjustable parameter u determines the 
strength of odd partial-waves between the baryons. 

(i) NN potential 

The parameters of the Minnesota potential were determined so as to reproduce the low-energy NN scattering data. 
This potential is adopted as the NN potential in what follows. One nice, noteworthy point is that the Minnesota 
potential reproduces reasonably well both the binding energies and sizes of few-nucleon systems such as 2 H, 3 H, 3 He 
and 4 He |13| . This is a necessary condition for the present purpose. Results of our calculation are insensitive to the 
choice of u because even partial- waves are predominant components, so u may be set to unity. The Coulomb potential 
is included in the present calculation. 

(ii) AJV potential 

Due to the experimental limitation of Ap scattering, the AJV interaction is rather poorly known. Though various 
versions of a model description for the AJV interaction can reproduce the existing Ap scattering data, the properties 
given by them are different from each other. For example, the 1 Sq and 3 Si phase shifts scatter depending on the 
model, though they are quite important for binding the A-hypernuclei. 

The relative importance of 1 Sq and 3 Si AN potentials may be estimated in a simple core nucleus+A model for the 
s-shell A hypernuclei The A separation energy, B\, may be estimated by 73a = —(A — 1)(Vr) — N S (V S ) — N t (V t ) 
when the kinetic energy contribution is neglected, where (V s ), for example, is the average of the AJV potential matrix 
elements in singlet states, and the number N s (N t ) of AJV pairs in the singlet (triplet) state is obtained by 

where J c is the spin of the core nucleus, Sa is the spin of A, and the total angular momentum J of the s-shell 
A-hypernuclei is J=J c +Sa- Table || lists the number of pairs obtained by assuming that the spin of the nuclear 
core remains unchanged by the addition of A. The fact that the ground state spin of A H is J = | but not J = | 
indicates that the 1 Sq potential is more attractive than the 3 Si potential. This conclusion is also consistent with the 
fact that the ground states of the A = A hypernuclei, A H and ^He, have J = instead of J=l, The energy splitting 
of the two states with different J values is expected to give us some information on the spin-dependent part of the 
AJV interaction. The above argument is qualitative. It is noted that in reality the dynamical effect of the AJV-EJV 
coupling modifies the contribution of the 1 S'o and 3 Si- 3 £>i AJV interaction components from the value of Table |, as 
shown in Ref. || for the A = 3, J=\ system. 

We adopt, for the AJV potential, the same form (|l|) as the NN potential. This ansatz is reasonable because of 
the similarity between NN and AJV, which is expected from the flavor SU(3) symmetry. In particular the 1 S'o AN 
channel is dominated by (A/i) = (22) SU(3) symmetry jf?]], which is exactly the symmetry of the 1 Sq NN channel. 
The CSB term of the AJV potential will be considered in subsect. 3.4. 

The potential depth and range were determined as follows. First we determined the potential parameters so as to 
reproduce the AJV S*-wave phase shifts predicted by the NF §| or the FSS (§ model (see Fig. @). Then the binding 
energies of the A = 3, 4 A-hypernuclei were calculated by using this phase equivalent potential. By taking a suitable 
combination of these two different potentials we attempted to reproduce the binding energy data. The iS-wave phase 
shifts calculated by the resultant potential was then fitted by a three-range Gaussian potential ([l]) . Of course the last 
procedure does not lead to a unique parameter set. In fact we could find a set of the parameters for a given value 
Vox of the repulsive part. Three sets of parameters determined in this way are listed in Table |l| 

The potentials of sets A and B, used in Rcfs. and have a moderate height of the central repulsion, 

while the repulsion of the set C potential is chosen to be exceedingly large. The potentials of quite different core 
heights are constructed in order to examine the extent to which the binding energy and the quark Pauli effect depend 
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on the potential shape, that is the dependence of the energy on off-shell properties of the potential. Though the 
potential shapes are different from each other as shown in Fig. |], they give practically the same low-energy scattering 
parameters, i.e., the phase shifts, the scattering lengths and the effective ranges. The S- wave phase shifts calculated 
by these effective potentials are compared in Fig. to those by the NF and FSS models. The 1 S and 3 Si phase shifts 
predicted by the effective potentials grow up to 32 and 19 degrees, respectively. This difference in the phase shifts 
is necessary to explain the energy splitting of the + and 1 + states of the A=4 systems. The AN potential used by 
Dalitz et al. || gives similar phase shifts between 1 Sq and 3 Si states. 

The parameter u of the AN potential is rather insensitive to the binding energy data of the A = 3 — 5 hypernuclei 
because the AN S-w&ves play a predominant role. This parameter was determined to be u = 1.5 so as to fit the 
forward and backward ratio of the Ap scattering data flT8|[l9"[ , as was done by Dalitz et al || . The sensitivity to u of 
the ratio is shown in Fig. |3|. 

The total cross section for Ap scattering is shown in Fig. ^. Our parameter set underestimates the experimental data 
in low-energy region. Some improvement is obtained by introducing the CSB term, as we will discuss in subsect. 3.4. 

Miyagawa et al. [|| have reported that the Nijmegen model NSC89 || reproduces the binding energy of ^H. 
According to their analysis, the So attraction of the AA^ interaction is quite important for reproducing the binding 
energy of ^H. The scattering lengths of the model NSC89 are a s = — 2.78 fm and at = —1.41 fm in a charge symmetric 
case. The corresponding lengths predicted by our effective potentials agree with those values within 15 %; a s — — 2.52 
fm and at = —1.20 fm. 

(hi) AA potential 

Since there are no data for AA scatterings, AA potentials are derived theoretically and the soundness of the potential 
may be tested against the binding energy data of the AA-hypernuclei. In the flavor SU(3) symmetry the dominant 
component of the AA channel is again (A/i) = (22) |l7j but its probability is reduced to about 70 %. Only three 
AA-hypernuclei have so far been confirmed experimentally. The binding energies of these nuclei seem to imply that 
the AA interaction is weakly attractive. The strength of the AA attraction in the AA-hypernuclei may be estimated 
by the Ai?AA value: 

AS AA ( A 1X) = Ba(aaX) - Ba( A ~a^) 

= B AA ( A ix)-2B A ( A ~lX). (3) 

The experimental AB AA values are AB AA ( A ^He) = 4.7 ± 0.6 MeV AB AA (l° A Be) = 4.3 ± 0.4 MeV §§ and 
ABaa(aa B ) = 4 - 8 ± °- 7 MeV |§- (Another interpretation of the last data implies A£ A a(aa b ) = ~ 49 ± °- 7 Me V.) 
We note in passing that the corresponding value for the nn attraction may be estimated as 

AB nn = B( 3 R) - 2B( 2 K) S 4.0MeV, 

which is in contradiction with the expectation that the nn interaction is more attractive than the AA interaction. 
This contradiction is probably because both 3 H and 2 H are rather weakly coupled systems. 

We use two types of AA potentials which predict weak attraction: One is the potential, denoted OBE-sim [ p5[ , 
which is phase equivalent to the Nijmegen model D (ND) potential ||, and another is the potential, denoted FSS-sim, 
which is phase equivalent to the quark cluster model FSS (26) . See Table [II] for the potential parameters of FSS-sim. 
The scattering lengths and effective ranges predicted by these AA potentials are 

a s = -2.80fm, r s = 2.81fm (OBE-sim) 
a s = -3.01fm, r s = 2.16fm (FSS-sim) 
a s = -3.01fm, r s = 2.14fm (FSS). 

Figure || compares the 1 Sq phase shifts calculated by the three models, and exhibits the potential shapes of the 
OBE-sim and FSS-sim models. 



B. Variational Trial Functions 

The Hamiltonian H of the system comprising nucleons and A particles is given by 

H = T + V = J2i^-T cm + V^ NN ^ + + + , (4) 

2—1 
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where is the Coulomb potential between the protons. The mass of N is taken as h 2 /uiN = 41.47 MeV-fm 2 

and the mass of A is set to be mjy/niN — 1.18826. In the present study, we calculate the binding energies of various 
systems up to A = 6 in a complete A-body treatment. The trial function for the eigenfunction of H is given by a 
combination of basis functions: 



K 



^JMTM T = CkPk, (5) 

fc=l 

where the basis function ip k is given by 

tpk = A{G(x, A k )xkjMVTM T }- (6) 

Here A is an antisymmetrizer acting on the identical baryons, x — (xi, • • • , xa-i) stands for a set of relative (e.g., 
Jacobi) coordinates, and XkSM s (Vtm t ) i s the spin (isospin) function of the system. The CG, G(x,Ak), is defined by 

A 

G(x,A k ) = exp j - y^QfcijQ-t - r 3 ) 2 | 

i<j 

1 A ~ X 

= exp { - - Ak v x i ■ x i)- ( 7 ) 

i,3'=l 

The (A— 1) x (A— 1) symmetric matrix A k contains A(A— l)/2 independent matrix elements which serve as nonlinear 
parameters to characterize the CG basis. Since 7*j — Vj can be expressed as a combination of x m , (r$ — rj) 2 becomes 
a combination of quadratic terms x m -x n , so the two expressions in Eq. (Q) are equivalent. The nonlinear parameters 
Akmn can be uniquely expressed in terms of a k ij and vice versa. The parameters (Xkij °f the former expression are 
more direct in representing the range of the correlation between the constituent particles |lj,^,^8) . 

The trial function (Q) has clearly the total orbital angular momentum L = 0. We assume that the hypernuclei 
treated in the present study are well described by only L = component. It is noted, however, that the partial waves 
corresponding to the coordinate Xi are not restricted to =0, but in general include higher orbital angular momenta 
[^9|,^0[. This is apparent because the basis function (Q) contains the cross terms Xi-Xj in the exponent of the CG and 
the expansion of those terms into polynomials contains high orbital angular momenta. 

The suffix fc in XkSM s is used to distinguish possible independent spin functions for a given S value |l3|]l4| ]. For 
example, in the case where four particles have the total spin 5'=0, two independent spin functions are possible. We 
took into account this possibility in the calculation. As for the isospin part we used an appropriately coupled function 
Vtm t for a given T value. 

The set of linear variational parameters c=(ci, • • • , ck) of Eq. (||) and the energy E are determined by Ritz varia- 
tional principle, which leads to the generalized algebraic eigenvalue problem 

He = EJVc, (8) 

where TL and Af are the matrices of the Hamiltonian and of the overlap 

Hij = (<pi\H\<pj), and N VJ = (<Pi\<Pj) = 1, • • • , K). (9) 

The matrix elements of the CGs can be calculated analytically for most of important operators. This plays a key 
role in obtaining a precise solution for few-body systems. Formulas for the necessary matrix elements are given in 



Refs. [14,271] and |28f|. See also the appendix where we show how to obtain the matrix element for many-particle 



density operators. 



C. Stochastic Variational Method 



The basis function tp k depends on A(A — l)/2 nonlinear parameters a ki - or A ki . and also a discrete parameter 
to specify the intermediate spin quantum numbers in the case when a number of independent spin functions, g, is 
g > 1 for a given S value. These parameters define the shape of the basis function and determine how well the 
variational function space contains the true eigenfunction. To find the best possible solution, one has to optimize 
the parameters. By assuming that we need a linear combination of K functions, we face an optimization problem of 
K(A(A — l)/2 + e) parameters, where e is when g = 1, or 1 otherwise. The number of parameters becomes very 
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large to get reasonably accurate solutions. For example, in the case of ^1^-6, we need K = 200 roughly and then end 
up with 3000 parameters. Moreover, Eq. (||) must be solved to calculate the energy, so that the optimization of the 
parameters is extremely hard in problems of the present type. 

Another serious problem of the minimization of a function is the omnipresence of local minima. A local minimum is 
the point where the function reaches a minimum in a finite interval of variables and the number of such minima tends 
to increase exponentially with the size of the problem. The conventional deterministic optimization algorithms (like 
the Powell or the conjugate gradient method) f3l| tend to converge to whichever local minimum they first encounter. 
To avoid such problems we use the SVM strategy. A key of the SVM is a one by one increase of the basis size by 
searching the best among many random trials for the basis function. 

The SVM strategy we have used here consists of two procedures, step-by-step and refinement |l^ , p^ , |3^| . The 
first stage is the following trial and error procedure to increase the basis dimension: Let us assume that the sets 
A\ , ■ ■ ■ , Ak~i and g\ , • • • , gk-i have already been selected, and the (k — l)-dimensional eigenvalue problem has already 
been solved. Here the index <?i specifies which spin function is chosen for ifi from among the g spin functions. The 
next step is the following: 

step-by-step 

1. Different sets (A\, • • • , A^) are generated randomly. 

2-1. For each set A\, the g eigenvalue problems of fc-dimension are solved with the g different spin functions and the 
corresponding energies (E\ , • • • , E^f) are determined. 

2-2. The parameter of the spin function that produces the lowest energy E l k among the set (Eg, ■ ■ ■ , E^ 9 ) is selected 
to be g\ for the set A l k . 

3. By repeating the above processes 2-1 and 2-2 from i=l to n, the energies (E^, • • • , E 7 k l ) are determined. 

4. The parameter set A™ and g™ that produces the lowest energy among the set (E\, • • • , E k ) is selected to be the 
fcth parameter set for (fk- 

5. Increase k to k + 1 and cycle the processes 1—5. 



The essential reason motivating this strategy is the need to sample different parameter sets as fast as possible. The 
advantage of this procedure is that it is not necessary to recompute the whole Hamiltonian matrix nor is it necessary 
to perform a new diagonalization at each time when a new parameter is generated This competitive selection 
substantially improves the convergence. 

Obviously the basis size cannot be increased forever. Moreover, when the A:th basis state is selected, the previous 
states are kept fixed, that is we tried to find the optimal state with respect to previously selected basis states, but 
actually some of the basis states selected earlier might not be so important anymore because the succeeding states 
took over their roles. So one may include a refining procedure where the previous states are probed again as described 
in the following: 

refinement 

1. Random parameter sets (A\, • • • , A£) are newly generated. 

2. The parameters of the fcth basis state are replaced by the new candidates and the energies (E\,, • • • , EZ-) are 
calculated. 

3. If the best of the new energies is better than the original one, then replace the old parameters with the new 
ones, otherwise keep the original ones. 

4. Cycle this procedure through the basis states from k = 1 to K . 



The selection of the spin function has to be performed in the step 2 in exactly the same way as in the step-by-step 
case. One may repeat the refining process till no further significant improvement is obtained. The refinement is useful 
for suppressing the basis size. 

In most of the calculations, we took n as n ~ 50 x (A(A — l)/2) and generated the matrix A^ by changing 



a kij (i < 3, i = 1; ' ' ' j A — 1) one by one randomly in the interval listed in Table HI, where, by expressing a as 

a = i (10) 



G 



d was chosen in the range of [0, d max }. Also listed in the table are values of g used in the present calculation. We 
show the convergence of the energy for ^iHe in Fig. |[ Both procedures of SVM are found to be very effective for 
obtaining solutions with high quality. 

A virial ratio 77 is often used to test the accuracy of the solution <3>. According to the virial theorem, the ratio 

A dV 

with W = Yr l -— (11) 

2—1 

must vanish for the eigenstate 'J of the Hamiltonian |14|] . We will calculate rj to check the accuracy of our variational 
solution. Note that a spherical symmetric potential V(\ri — rj\) is changed into rV'(r) with r = |r< — Tj\ by the 
operation of J2i=i r i' W~- 



(if>\W\if>) 
2(*|T|*> 



III. RESULTS 



First we show in Table IV the root-mean square (rms) distances of the particles in the A- and AA-hypernuclei 
together with the typical dimension size K and the virial ratio rj. Here the rms distances are defined as follows: 



(<■*) 



3<* 



i=l 



), (r2.)=(*|(r,-r,) 2 |*), 



(12) 



where xa is the center-of-mass coordinate of the system. We see that both and \fR, amon g others, have very 
large rms values, reflecting the very small binding energies, as shown below. The rms distance of NN decreases by 
adding more A particles: By adding A to 2 H, the {r 2 NN ) value changes from 3.9 to 3.6 ( 3 H) and 3.3 fm ( A ^H), 
respectively. A similar but slightly mild change is seen in 3 H and 3 He as well. In the case of 4 He, however, we see 
that the rms distance of TV TV hardly changes by the addition of A. Generally speaking, the basis size K needed for 
an accurate solution increases with an increasing number of particles and also with an increasing spatial extension 
of the system. Comparing the solutions for the different potentials, the potential with an exceedingly large repulsion 
like the set C or the OBE-sim potential makes an accurate solution more challenging. It is worthwhile noting that the 
increase of the basis size is rather gentle on these conditions, so that the solutions have been obtained in the order of 
rj < 10~ 4 in most cases. 



A. A-Hypernuclei 

Table ^ shows the results of calculation for the A-hypernuclei. Since the parameters in Table |J are determined 
to fit the binding energies of the A— 3, 4 A-hypernuclei, it is apparent that the energies of those nuclei are in good 
agreement with the experimental data including the excited states of the A=A system. As noted in subsect.2.1, we 
expect that our effective potential has very reasonable relative strength between 1 5o and 3 <Si states. 

Though the potentials of set A and set C are quite different in the shape, they are constructed to be phase equivalent, 
and the binding energies predicted by these potentials are almost the same. It is an agreeable feature that the phase 
equivalent potentials predict the same energies for all the nuclei. This is a strong point of the present calculation 
with high precision. Of course the solution for the case of set C potential requires a larger basis dimension than the 
one for the case of set A potential, so that the variational solution must be obtained to an accuracy of a few tens of 
keV regardless of the characteristics of potentials. We note that the above nice feature is not preserved in a naive 
calculation such as a frozen core nucleus+A model and that you may draw an erroneous conclusion based on such 
calculations. For example, if we describe ^He as an a + A model where the a-particle wave function is given by a fixed 
(0s) 4 Slater determinant with the size parameter of 1.39 fm, we obtain quite different B\(\}ie) values depending on 
the core height of the A7V potential: 3.95 and 2.29 MeV for set A and set B, respectively (no anomaly in _BA(iHe)!). 
Set C potential does not even bind ^He. (Here B\(jjie) is defined by the difference of the energy of the model 
a-particle (—23.94 MeV) and the total energy of the a + A model calculation.) 

As seen in Fig. [l], our effective potential is more attractive in 1 Sq than in 3 Si. We asked a question of whether this 
potential is strong enough to bind an excited state of ^H with J 7 * = | + . The calculation showed that the energy of 
the system approaches the deuteron energy, that is, the potential accommodates no such a bound state. 

It has been known for a long time that the CSB effect of the AiV interaction is manifestly exhibited in the Ba values 
of the A=4 hypernuclei. The data suggest that the Ap interaction is stronger than the An interaction. Though the 
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AN-T,N coupling or A-S° mixing is expected to play the significant role in the CSB, no conclusive understanding of 
the mechanism of CSB has been obtained yet. Both potentials of sets A and C are charge symmetric. An extension 
to a calculation including the CSB effect will be discussed later. 

The result for A He confirms the long standing problem, i.e., the potential fitted to the A— 3, 4 data overbinds A He 
by about 2 MeV. We will come back to this anomaly in sect. 4, where the effects of the quark Pauli principle are 
discussed. 



B. AA-Hypernuclei 



Table VI shows the results of the AA-hypernuclei. The energies calculated with different AA~ potentials are almost 
the same. We have examined if a lightest system of A — 3 forms a bound state, but have not found a stable AAA 
system with the AA potential either OBE-sim or FSS-sim. We have then asked a question of whether or not a 
bound state is formed for A > 4 system. We have calculated the energies for various spin S and isospin T states: 
(S, T) = (0, 0), (0, 1), (1,0) and (1, 1) for A ^H, and (S, T) = (0, 1) and (1, 1) for A ^He or AA m. We find that a lightest 
bound state of AA-hypernuclei is only AA H with (S,T) = (1,0) in consistency with the conclusion of Ref. [g5|. See 
Table [VJ. The binding energy is so small that its existence may be strongly subjected to the strength of the AA 
potential. This state has a structure of a deuteron plus AA that is coupled to S—0. Results for A— 5 AA-hypernuclei 
are also listed in Table [v| It is very likely that both AA H and AA He are bound. Experimental search for these 
AA-hypernuclei and measurements of their binding energies will provide us with valuable information on the AA 
interaction. 

The main component of AA is the 1 Sq state for the A=4— 6 AA-hypernuclei, so the total angular momentum J of 
each hypernucleus is equal to that of the core nucleus. The J and T values are (J,T) = (|, |) for AA H and AA He, and 
(0,0) for AA He. The contribution of the 3 Si and the 1 S'o components of the AN potential to the potential energy is 
determined by N t and N s values. Using an arithmetic similar to Eq. (|^) we obtain N t —^(A — 2) and N S =^(A — 2). 
Their ratio is Nt : N s = 3 : 1 for all these AA-hypernuclei, as in the case of A He. 

The energy of BA\( A ^He) is again overbound compared to the empirical value, just as in the case of £?a(aHc) 
shown in Table [y| The calculated value of A£?AA( AA He) is, however, in reasonable agreement with the experimental 
value of 4.7 ± 0.6 MeV: AB AA = 4.3 MeV for OBE-sim AA potential and 5.2 MeV for FSS-sim AA potential. We 
thus expect that the strength of the AA potentials employed in the present study is nearly right. Rather we guess 
that the same anomaly known in A He appears in AA He as well. 



C. Density Distributions and Correlation Functions 



Obviously the density distribution and the correlation function between the particles provide us with more profound 
information on the structure of the system than just the average distances of the particles. It is also of particular 
interest to know the extent to which the core nucleus is distorted by adding the A particles. For this purpose we 
define the distribution function of A?" or A by 

p(r) = <*| < y(r i -fl c -r)|*), (13) 

where r j denotes the position vector of N or A, and R c is the center-of-mass coordinate of the core nucleus. Therefore 
r is a position vector from the core nucleus. The density defined above is different from the conventional one which 
is calculated as a function of a position vector measured from the center-of-mass coordinate of the system. The N or 
A density defined here seems more direct to represent the change with the addition of A. See the appendix for the 
method of calculation of the density and the correlation function etc. 

Figures @ — || show the density distributions of N and A when A particles are added to the core nucleus of 2 H, 
3 H or 4 He, respectively. The set A AN and OBE-sim AA potentials are used but other choices of the potentials 
do not change the result very much. As seen from these figures, the density distributions of N and A are quite 
different. The A distribution reaches further more outwardly than the N distribution. The systems of A H and AA H 
are extremely spread out because the deuteron is weakly bound and in addition these hypernuclei have very small 
separation energies. The maximum probability of finding A occurs at or near the point where the nucleon probability 
falls to half of its maximum probability. Figure || shows that the A distribution is considerably different between the 
ground and excited states of A H and that the core distortion is larger in the ground state than in the excited state. 
We see that the addition of A particles makes the core nucleus shrink and accordingly the A density move to the inner 
region. 

The function C(r) defined by 
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C(r) = {y\5(r i -r j -r)\V) 



(14) 



gives the information on the correlation of the constituent particles. The correlation function of AA is calculated 
with the set A AA?~ and the OBE-sim AA potentials and compared in Fig. [jTJ for the AA-hypernuclei, aa-^' aa^ 
and AA-He. It is noted that the correlation function has vanishingly small amplitudes at short distances, as it should 
because the OBE-sim AA potential has a strong central repulsion. Since the CGs have large amplitudes at short 
distances of the particles, this vanishing of the correlation function indicates that the CG basis is sufficiently flexible 
to reproduce such characteristics and that the SVM has been practical in selecting suitable basis sets. The behavior 
of the correlation function reflects the ABaa value defined by Eq. (||): The calculated A_Baa values are 0.06, 1.2 and 
4.3 MeV for aaH' aaH and A ^He, respectively. The wider the spatial extension of the correlation function, the weaker 
the AA interaction. We see that the two A particles, particularly in aa-H' spend almost all the time outside their 
interaction range. The existence of AA H crucially depends on the strength of the AA potential. The peak position of 
the correlation function moves depending on the size of the underlying core nucleus. 

To visualize the structure of the system, it is useful to calculate the two-particle (particularly AA of the AA- 
hypernuclei) distribution function given by 

D(r,r') = (y\6(n -R c -r)S(r 3 - i? c - r')| *)• (15) 

This function depends on the three quantities, r, r' and 9, in the case of L = wave function, where 9 is the 
angle between r and r'\ r r' = rr' cos 9. In Fig. O, the two-A distribution function D(r,r') weighted by r 2 r' 2 sm6 
is displayed for aaH- To draw this function we first searched for a point where the probability, r 2 r' 2 sm9 D(r, r'), 
reaches a maximum. The maximum for aaH appears at the point of r = r'=1.8 fm and 6*=47°. Drawn in Fig. |ll| is 
the distribution of A when the other A was placed at that point. Though the point determined above may suggest a 
configuration of an isosceles triangle with side lengths of 1.8, 1.8 and 1.4 fm, the distribution is actually spread out 
in a wide region, so that the geometry is to be understood rather loose. For example, the rms distance of AA is as 
large as 3.5 fm, much larger than 1.4 fm. 

The core nucleus can be distorted by the addition of A particles, that is, the nucleus never remains in its ground 
state in the hypernucleus but may have mixing in of other states. To calculate the probability that the core nucleus 
remains in its ground state, we define the spectroscopic factor S by 

S = E / dr{g McMs {r)}\ (16) 

Me, M s 

where the spectroscopic amplitude gu c M s (?") for the A-hypernuclei is given by 

<7M C M» = (^J aMc TM T { A ^)XlM s {^{xA-l - f)|tf JMTM T {i*)) ■ (17) 

Here xa-i = va — R c is the distance vector of A from the center-of-mass of the core nucleus A_1 X, and 
"i 1 j c m c tm t { A ~ 1 X) is the ground state wave function of the core nucleus. If no distortion of the core nucleus is 
present in the hypernucleus, Y] M m {9m c M, ( r )} 2 reduces to the density distribution of A that is defined by Eq. (|l3|). 

The spectroscopic factor for the AA-hypernucleus is calculated in a similar way. The spectroscopic amplitude may 
be defined by 

g Mc (r,r') = (^j c m c tm t ( A ~ 2 X) X oo(AA)(5(a^-2 - r)6(x A -i ~ r')\^ JM TM T (aa X)), (18) 

where xa-i = r A —i — Rc and xa-i = t a — Rc+a with R c +a denoting the center-of-mass of the core nucleus and 
the A. Note that the two A particles were assumed to be in spin singlet. The spectroscopic factor S is then obtained 
by {ffMefcj*')} 2 summed over M c and integrated over r and r'. 



Table |VII| lists the spectroscopic factors calculated for A- and AA-hypernuclei. We confirm, in conformity with 
the results of Figs. [?| — that adding more A particles results in smaller S values, that is, it leads to the larger 
distortion of the core nucleus, as expected. There is no experimental information on S values but just one theoretical 
value is available for A H. According to Ref. || the value is 0.987, which is in very reasonable agreement with the 
present estimate (0.991) in spite of the use of quite different potentials. The S values of A He and AA He are fairly 
large compared to those of the other nuclei, which may indicate that the 4 He core is considerably stable. As already 
pointed out, however, these results do not necessarily substantiate that the frozen core model is a perfect model even 
for the light s-shell hypernuclei. A small distortion becomes sometimes crucial for the accurate evaluation of the 
binding energies. 
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D. Charge Symmetry Breaking 



The experimental data of the A=4 system imply the CSB of the AN interaction. The CSB is observed for both the 
ground and excited states: 

AB A = B A ( A Ke) - Ba(aH) = 0.35 ± 0.07MeV, 

AB A = B A ( A He*) - B A ( 4 A R*) = 0.24 ± 0.15MeV. (19) 

Bodmer and Usmani J3(| investigated the CSB phenomenologically. They suggest that the CSB interaction is effec- 
tively spin independent, that is, it has no <t a -<Tn dependence. In their argument, however, the ground and excited 
states were not treated simultaneously. They modified the strength of the AiV attraction with respect to the ground 
or excited state, in spite of the fact that the strength of <t a ■ cr^v interaction is different from each other in the 
case of J 71 = + and 1 + states. Dalitz et al. || studied the CSB effect for both cases of the spin-dependent and 
spin-independent CSB potentials. However, they focused only on the CSB effect of the ground states of \H and 
A He. Since the accuracy of the present calculation is high enough to discuss the CSB, we attempted to determine 
phenomenologically a AN interaction which includes the CSB effect. 

We introduce the CSB term in both components of the triplet and singlet AN potentials (|j) as follows 

V ot -> V ot - v t T 3 , V Qs -> Vos - v s t 3i (20) 

where t 3 is the third component of the isospin matrices: 

r 3 \n) = +\n), r 3 \p) = -\p). (21) 

The modification ( p0[ ) of the strength implies the CSB potential of the form 

Vcsb = {\( 3v t + v s )t 3 + \{v t - v s ){ct a -(t n )t^ e~ Kr2 . (22) 

(Here K t and k s are assumed to be the same to simplify the expression but no such an assumption was made in the 
calculation.) If v t and v s have the same sign, the spin dependence of the CSB potential becomes weak. 

Starting from the set A AA^ potential in Table 0, we determined those Vt and v s values of the CSB potential 
which reproduce the data as suggested in Eq. (|l9|). For this purpose we examined the energy dependence of the A=A 
hypernuclei on the set of vt and v s values around (vt,v s )=(Q, 0). After repeating calculations for different sets of 
parameters, we found a set of values which reproduces both AB A and AB A very well: 

v t = 3.30 MeV and v s = 2.65 MeV. (23) 

Both vt and v s values are positive in sign, so that the CSB potential determined above is weakly spin-dependent 
in accordance with the conclusion of Ref. With this CSB potential the Ap phase shifts in 1 S'o and 3 Si states 

increase to 34 and 21 degrees, respectively. 

As listed in Table |y|, the calculated binding energies with this CSB potential are in good agreement with experi- 



mental data. The change of B AA with the CSB potential is listed in Table VI. The CSB does not alter the prediction 



for possible existence of and A ^He. 

The CSB potential produces a very minor change in energy for those hypernuclei which include the same number 
of protons and neutrons such as A H, AA H, ^He and AAHe. 



The scattering parameters calculated with the CSB potential (|23| ) are compared in Table VIII to predictions by 
different models. Some models by e.g. Nijmegen group P,p|, ^0| , p7| give a spin-dependent CSB interaction and predict 
the strengths of the AJV attraction in the order of 

|aM > \a^\ > |af p | > \4 n \. (24) 

On the contrary, our potential predicts the opposite order in the So part: As expected, the Ap interaction becomes 
more attractive than the An interaction. In the case of 3 5i both the Nijmegen and our phenomcnological potentials 
predict that the Ap interaction is more attractive than the An interaction. It is noted that the introduction of the 
CSB determined in this study improves a fit to the Ap total cross section as displayed in Fig. 0. 
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IV. QUARK PAULI EFFECT 



A. Pauli-forbidden Configuration 



As shown in the previous section, the calculated A separation energy of jjHe is overbound by about 2 MeV, 
compared to experiment. The AA separation energy of ^^He is found to be also largely overbound. To resolve these 
discrepancies is one of the most challenging problems in s-shell hypernuclei. Shinmura et al. |Q argued that the 
noncentral components of the AN and AW potentials plays an important role. The effect of the A-S conversion 
may become important, as shown by Miyagawa et al. ||, through the 3 S- 3 D tensor coupling in the AN-T,N channel. 
An ab initio calculation using realistic potentials has not been performed until now because of the complexity of the 
potentials. Even if it has been done in some way, there is so large uncertainty in the realistic potentials that such a 
calculation may not give a perfect answer. The premise we are taking in this study is that the effects of the noncentral 
forces and the A-S conversion etc. are renormalized in the effective potential which is determined so as to fit the 
binding energies of the A— 3, 4 A-hypernuclei. We thus pay attention to other mechanism, a quark Pauli effect, for 
examining this problem in this section. The quark Pauli effect was considered in Ref. |38| on the basis of the quark 
shell model where the quarks were allowed to move in the whole hypernucleus. In contrast to this approach, we treat 
the Pauli effect by confining three quarks inside the baryons. Parts of the results along this line have been reported 
in Refs. Q and Q. 

So far we have treated the A particle as distinguishable from the nucleon N. Therefore all the baryons in ^He or 
A ^He, for example, may occupy a small region at the same time if that configuration is favorable to the energy gain. 
In the quark model of baryons, however, A and N share quarks of the same flavor, and then in a A-hypernuclei with 
A > 5 there may be an increasing chance that more than six u or d quarks occupy the same orbit as the baryons 
come close to each other. Of course this chance is forbidden from the quark Pauli principle, so that it should be 
excluded in a calculation 

011 • This type of quark Pauli principle should act even in normal nuclei but its effect is 
taken into account by the Pauli principle for the N. There is no way, however, to take into account the quark Pauli 
effect in those calculations which do not consider the quark substructure of the baryon. The aim here is to estimate 
how much the A separation energy changes by the mechanism arising from the quark Pauli effect. We perform this 
estimation by using only those wave functions which are still expressed in terms of the baryon coordinates. The effect 
to be discussed here is concerned with a specific forbidden configuration, which is a special effect of a more general 
antisymmetrization effect of the quarks |39[| . 

A specific quark Pauli-forbidden configuration comprising five baryons was considered in Refs. [ pd| and [ O ]: 

m*) = (^)" ex p {-i X> - ^) 2 } > ( 25 ) 

where R5 is the center-of-mass coordinate of the five baryons. This configuration was derived by assuming that a 
baryon is a three-quark system and that the orbital motion of each quark is described by a (Os) harmonic-oscillator 
function, exp{— r 2 /2b 2 }. The parameter b determines the size of the baryon. If the simple (Os) 3 wave function ought 
to reproduce the charge radius of proton (0.86 fm) jh|, then b becomes 0.86 fm. A more reasonable value of b would 
be slightly smaller than 0.86 fm if other effects such as meson clouds etc. are considered in the charge radius. We will 
examine the binding energy change in the range of 6=0.6—0.86 fm. 

Now the trial function must be subjected to remain in the subspace which is orthogonal to the forbidden state, that 
is, (V'fI*) = 0. We thus modify Eq. (§) by 

^ JMTM T = Cfc I Ifk — £r a <r tt |v> fc > , (26) 

fc=l \ a / 

where T a is an orthonormal set which spans the Pauli-forbidden space. In the case of A=5, T a can simply be 
represented by AlipFXV}- F° r > 6, however, the construction of T a is never trivial. The system A ^He, for example, 
has three types of forbidden states ( |25| ) comprising ppnnA, ppnAA, or pnnAA baryons. A method proposed in Ref. [12 ] 
enables us to eliminate practically the Pauli-forbidden components. The reader is referred to Ref. |T|] for its detail. 
The basis functions ipt were newly selected by the SVM. The accuracy of the solution is an essential ingredient in 
this calculation because the solution must be accurate enough to be sensitive to small components produced by the 
quark Pauli principle. 
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B. Hypernuclei with A=5 and 6 



Table IX shows the dependence of the A separation energy on the baryon size for the A=5, 6 hypernuclei. The case 



of 6=0 is no quark Pauli calculation. The energy change is insignificant for b < 0.6 fm. As expected, a significant 

A He and AA I 

AA H and AA He. The energy of AA He is apparently overbound in no quark Pauli calculation. The elimination of the 
forbidden states produces a surprisingly large reduction in the binding energy, the order of half a MeV even at &=0.6 

5t 
AA 1 



binding energy reduction is found to occur at b = 0.86 fm in A He and AA He, while the reduction is much less in 

■n i 

at 

fm, which is very favorable to the anomaly problem. These results can be understood from the fact that AA H and 

8 
AA 



He are spatially more extended than A He and AA He. (See Table and Figs. [?] — Since the ^-dependence of the 
quark Pauli effect is very moderate in AA H and AA He, a measurement of their binding energies will be quite useful 
for the determination of the AA attraction. 



As the forbidden state (|25|) is a spatially compact configuration, one may conceive that the energy change be 
dependent on the characteristics of potentials, especially the height of the central repulsion. Table |x| compares B\ 
and B\a values for different sets of the potentials. The variation of the A separation energy for the different potentials 
is very moderate, so we may conclude that the results of these tables indicate the general trend of the quark Pauli 
effect in light A- and AA-hypernuclei. 

To explain the mild dependence on the potentials we derive an approximate expression for the energy change as 
follows. Let and Eq denote the ground state wave function and the energy of no quark Pauli calculation. When 
the quark Pauli effect is taken into account, the ground state wave function changes and, since the change is expected 
to be small, the wave function may be approximated as 



The energy change is given by 



f= . - (Vo-siPf) with £ = (Vf|#o>- (27) 

yl - e z 



E-E Q = (* -Eon e 2 {^ F \H\^ F ), (28) 



where we used H^q = Eq^q and assumed that e is small and (^fI-HI^f) 3> Eq. The energy change is approximately 
estimated by the product of the probability, e 2 , of finding the forbidden component in and the energy expectation 



value of the Pauli- forbidden state, (V'fI^IV'f)- We show in Table XI e , (ipF\H\ife), and the energy change calculated 
with the different potentials. In the case of A He, for example, the set A potential with soft central repulsion gives 
e 2 = 0.00406, (ip F \H\ip F ) = 513 MeV, while the set C potential with very large repulsion gives e 2 = 0.00331, 
(ip F \H\ip F ) = 969 MeV. We see that these two factors really depend on the potential but their product is rather 
weakly dependent on the potential, and close to the binding energy change SE. This feature is a basic reason for 
the mild dependence of the energy change on the potential. A similar result is seen in the case of AA H, where the 
OBE-sim or FSS-sim AA potential produces almost the same energy change though the core height is quite different 
between the two potentials. It is also interesting to recognize the difference of the energy change between A He and 
AA H. Comparing the results calculated with set A and FSS-sim potentials we understand that, though (ip F \H \ip F ) is 
almost the same in both cases, e 2 is significantly different in both nuclei, which explains the larger quark Pauli effect 
in A He than in AA H. 



V. SUMMARY 



We have determined phenomenologically the AN central potential so as to reproduce the binding energies of 
A = 3, 4 A-hypernuclei on the premise that various effects such as the noncentral forces and the A-E conversion etc. 
are renormalized in the effective potential. The phase shifts predicted by this potential reaches 32 and 19 degrees at 
maximum in 1 Sq and 3 5i states, and the scattering lengths and effective ranges are a s = —2.52 fm, a t — —1.20 fm, 
r s = 3.08 fm and r t = 4.26 fm, respectively. 

Bound state solutions for the Schrodinger equation have been obtained by the stochastic variational method with 
the correlated Gaussian basis. We have performed a high precision calculation to determine the binding energies of 
the hypernuclei with A < 6. We have properly taken into account both the short-range correlation of the particles 
and the asymptotic behavior of the wave function in setting up the variational trial functions. We have employed the 
three phase-equivalent AN potentials which have different core heights and confirmed that they all give almost the 
same A separation energies for the s-shell A-hypernuclei ( A H, A H, A He, A H*, A He*, A He). 

The charge symmetry breaking part of the AN potential has been determined so as to reproduce both AB A and 
A£> A values for the A = 4 A-hypernuclei. We find out that the CSB is weakly spin-dependent in accordance with 
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the conclusion of Ref. |56|. The Ap potential is more attractive than the An one in both 1 5o and 3 Si states. The 
Ap phase shifts increase by 2 degrees with the inclusion of the CSB effect. This phenomenology disagrees with the 
prediction of some OBE models where the CSB is spin-dependent and the An potential is more attractive than the 
Ap potential in state. 

For the study of the AA-hypernuclei we have used the two different AA potentials which reproduce the Ai?AA( A ^He) 
value reasonably well. The present calculation has shown that we can expect hitherto undiscovered particle-stable 
AA-hypernuclei, aa^, aa-H and AA-He, with the Baa values of 0.4, 5.5 and 6.3 MeV, respectively. No other bound 
AA-hypernucleus has been obtained for A < 5. Since the binding energy of A a^ ^ s verv small, its existence will 
crucially depend on the attraction of the AA interaction, while both of aa^ an< ^ AA^ e ar e more tightly bound, so 
that their existence will be very probable. Experimental confirmation of these AA-hypernuclei will provide us with 
valuable information on the AA interaction. 

The structure change produced by adding A particles has been investigated by calculating the density distribution, 
two-particle correlation function and spectroscopic factor etc. The distribution of A is in general much broader than 
that of N. For instance, the root-mean-square distance of the two A particles in AA ^H is larger than 8 fm and they 
spend almost all the time outside the range of the AA interaction. 

The present calculation has confirmed the well-known problem that the binding energies of both a He and ^^He are 
predicted to be overbound by 2 — 3 MeV. We have focused our attention on the effect of the quark substructure of N 
and A to see if the binding energy is reduced by the quark Pauli principle. The effect is negligible if the baryon size 
in which the quarks are confined is smaller than 0.6 fm, but becomes appreciable, particularly in aa-H^ ^ the size is 
taken to be as large as 0.7 fm. To resolve this over binding problem still remains very challenging and important in 
few-body hypernuclear systems. 



APPENDIX A: MATRIX ELEMENTS OF THE CORRELATED GAUSSIANS 



Most of the formulas for the matrix elements of the CG are given elsewhere |14| , p7| -p9|. The reader is referred 
to these references for the detail, so here we recapitulate the basic formulas and the calculation of the correlation 
function. We consider the system comprising N particles and denote the set of relative coordinates as (xi, • • • , ccjv-i)- 

The CG has two noteworthy advantages in the calculation of matrix elements in addition to its flexibility to represent 
various shapes of many-variable functions. First the CG does not change its form under the permutation of particle 
indices. Since the permutation P induces a linear transformation of the coordinates as 

N-l 

P Xl = T H x i {i = N-l), (Al) 

we notice that the CG is changed by P into 

PG(x, A) — G(x, TAT) with T^ — T. }i . (A2) 

The action of P is thus very simple: The CG remains a new CG with the matrix A being replaced with a new one 
TAT. Thus the antisymmetry requirement for the basis function involving the CG is relatively easily fulfilled by just 
knowing the matrix T. We do not need to calculate other matrix elements than those between the CGs. Secondly 
most of the matrix elements for the CG can be obtained analytically, and this makes it possible to obtain a precise 
solution with the CG basis. 

The overlap of the CG is given by 

/'(2tt) n - 1 \ % 

(G(x,A)\G(x,A')) = [ i-i with B = A + A'. (A3) 

\ det B J 

By expressing the kinetic energy operator as 

— T rm = — > A»„-7r,' • 7T.; with 7T, = — , (A4) 

2m, cm 2 ^ 3 1 3 3 dxi' y ' 

we obtain the matrix element for the kinetic energy operator: 
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N 2 

(G(x,A)\^g--T cm \G(x,A')) 

1=1 

= hi 2 Tr {A' B- 1 AA) (G(x,A)\G(x,A')). (A5) 
To calculate the matrix element for the two-body interaction, it is convenient to express the potential in the form 

V(n -rj) = J drV{r)S(n - rj - r) (A6) 
and note that the relative distance vector of the particles can be written in terms of the coordinates x as 

JV-l 



^2 ^kXk = ojx. (A7) 



fe=i 

Then the matrix element of the potential is given by 

(G(x,A)|V(r < -r i )|G( a! ,A')) 
= / drV(r)(G(x,A)\S(uix-r)\G(x,A')) 



where 



{G(x,A)\G(x,A')) J drV(r)exp(-±c 



. (A8) 



- = uB^uj. (A9) 
c 

To prove this formula, one makes use of the Fourier integral for the (5-function and the many- variable Gauss integration. 
The advantages of this method are that, since the dependence on the form of the potential appears only through the 
last integration, different types of potentials can be treated on equal footing and also that the correlation function 
(O) can already be calculated in the above procedure. By choosing, e.g., V(r) = r" we can get the matrix element 
oF|rj — Vj\ n very easily. 

The calculation of the density distribution of the particle as defined by Eq. (|l3| ) poses no problem if one notices that 
Ti — R c is expressed in terms of a lin ear combination of (xi, • • • , x^—i). One only needs to redefine uj appropriately 
and then use the formula of Eq. (A8). 

By generalizing the method presented above for the matrix element of the (5-function, we can calculate the two- 
particle distribution function (|15| ) or even more-particle distribution function. All of these functions can be obtained 
through the following formula: 



(G(x,A)\l[6(wWx- Pi )\G{x,A')) 



i=l 

detC \ ' / 1 



(27T) n 



) CX P \-\ J2 C v Pi 'Pi 



i,3=l 

x (G(x,A)\G{x,A')}, (A10) 

where C is an n x n matrix defined by 

(C-% =wV)B- 1 w( j l (All) 
By integrating the above equation over p x , ■ ■ ■ , p n , one recovers the overlap (|A3|) of the CG. 
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FIG. 1. The 1 So (left) and 3 Si (right) phase shifts of AN scattering as a function of A momentum pa- Solid lines are obtained 
by the potentials of sets A, B and C, + by FSS, and x are by the Nijmegen model F, respectively. 

FIG. 2. The AA potentials of sets A and C as a function of AT-A distance r. 
FIG. 3. Forward-backward ratio of Ap scattering as a function of A momentum pa- Experimental data are taken from 

Refs. Ill and ff|. 

FIG. 4. TotaTcross section of Ap scattering as a function of A momentum pa- The dotted line is obtained with the potentials 
of sets A^B and C, while the solid line is obtained with set A potential including CSB. Experimental data are taken from 
Refs. (SJfQ. 

FIG. 5. The 1 5'o phase shifts (left) calculated by the AA potentials displayed on the right panel. 
FIG. 6. The energy convergence of Aj\He as a function of the basis dimension K. Set A AA" and FSS-sim AA potentials are 
used. The CSB effect as we discuss in subsect. 3.4 is included. The upper panel is the step-by-step selection up to K=2G0, 
and the lower panel shows the refinement at AT=200, which is followed by the step-by-step basis extension to A"=300 and the 
second refinement at A"=300. 

FIG. 7. The density distributions of N and A for 2 H, A H and aaH as a function of r, the distance from the center-of-mass 
of 2 H. 

FIG. 8. The density distributions of N and A for 3 H, A H, A H* and AA H as a function of r, the distance from the center-of-mass 
of 3 H. 

FIG. 9. The density distributions of N and A for 4 He, A He and AA He as a function of r, the distance from the center-of-mass 
of 4 He. 

FIG. 10. The correlation function of two A particles for the AA-hypernuclei, as a function of r, the distance between A-A. 
FIG. 11. The two-A distribution function r 2 r' 2 s'm8D(r, r') (upper) of AA H as a function of r and 0, where r' = 1.8 fm and 
(x,y) = (r cos 6, r sin 9). The core nucleus 3 H is located at the origin and one of the two A particles is placed at the point 
denoted X. The lower panel plots the contour line of the distribution function. Set A AA and OBE-sim AA potentials are used. 
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TABLE I. Contributions of 1 So and Si AN potentials to the A separation energy in a core-nucleus+A model. J is the total 
angular momentum and J c is the angular momentum of the core nucleus. 







J 




Jc N s N t 






i 

2 




1 2 I 

2 2 




3 
2 




1 2 


iH, iHe 









1 3 3 

2 2 2 




1 




115 
2 2 2 


iHe 




1 

2 




1 3 


TABLE II. Parameters of NN, AN and AA potentials. Three sets for the AN potential are charge symmetric, giving the 
same 1 So and 3 Si phase shifts in low energy region (pa < 600 MeV/c): The scattering length and effective range are a s — —2.52 
fm, r s — 3.08 fm, and at = —1.20 fm, rt = 4.26 fm. To include the charge symmetry breaking for set A AN potential, increase 
(decrease) both Vot and Vo s by the values denoted CSB for Ap (An). 






V R 

(MeV) 


Vot 
(MeV) 


Vos KR Kt Ks 

(MeV) (fm- 2 ) (fm" 2 ) (fm" 2 ) 


NN 


Ref. [16| 


200.0 


178.0 


91.85 1.487 0.639 0.465 


AN 


Set A 
CSB{t> 


200.0 


109.8 
{±3.3} 


121.3 1.638 0.7864 0.7513 
{±2.7} 




Set B 


600.0 


52.61 


66.22 5.824 0.6582 0.6460 




Set C 


5000.0 


47.87 


61.66 18.04 0.6399 0.6325 


AA 


FSS-sim 


200.0 




130.8 2.776 1.062 


TABLE III. Values of d 
possible spin functions are 


max, given in units of fm, and the number of spin functions g used in the present calculation. All 
taken into account except for AAHe, which is restricted to [[[S p , S p ]o, [S n , S n ]o]o, [5a, SaJoJoo- 


dmax Q 


2 H 
3 H 
3 He 
4 He 

in 4 

A H, a He 
A H*, A He* 
IHe 
A |He 

aaH, AAHe 
AAHe 








20 1 
15 2 
15 2 
12 2 
45 2 
18 2 
24 3 
15 5 
55 3 
18 5 
15 1 
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TABLE IV. The root-mean-square distances of A- and AA-hypernuclei. K is the basis dimension and 77 is the virial ratio 
defined by Eq. (|ll|). Set A AN and OBE-sim AA potentials are used. 





v( r > 


/ / 2 \ 
V( r NN) 


/ / 2 \ 


ill \ 


K 


V 




1.95 


3.90 






30 


2 x 10~ li5 


H 


1.71 


2.95 






30 


1 x 10~ 5 


3tt 

J He 


1.74 


3.01 






30 
60 


1 x 10~ 4 
6 x 10~ 6 


4u„ 

lie 


1 A 1 

1.41 


on 
2.30 






60 


1 X 10 


a tt 
aH 


4.9 


3.6 


10. 




30 


3 x 10~ s 


4 TT 


2.0 


2.7 


3.8 




100 


6 x 10~ 5 


aH* 


2.4 


2.8 


4.7 




100 
200 


2 x 10~ 4 
2 x 10~ 5 


A TT 

iHe 


2.0 


2.8 


3.8 




100 


1 x 10~ 4 


iHe* 


2.4 


2.9 


4.8 




200 


2 x 10~ 5 


iHe 


1.6 


2.3 


3.0 




200 


6 x 10 -5 


a|H 


4.2 


3.3 


7.1 


8.4 


200 


5 x 10 -4 


aaH 


2.0 


2.6 


3.2 


3.5 


200 


2 x 10 -4 


AAHe 


2.0 


2.6 


3.3 


3.5 


200 


2 x 10~ 4 


AAHe 


1.6 


2.2 


2.6 


2.8 


100 
200 


2 x 10~ 4 
8 x 10~ 5 



TABLE V. A separation energies and excitation energies, given in units of MeV, of A = 3 — 5 A-hypernuclei. a) Ref. |33j], b) 
Ref. H 





#a(aH) 


Sa( a H) 


B A (\Ke) 


S A (iHe) 


£ x (iH*) 


-Ex (IHe*) 


Set A 


0.18 


2.24 


2.20 


4.98 


1.14 


1.13 


with CSB 


0.18 


2.05 


2.40 


4.98 


1.08 


1.19 


Set C 


0.17 

n 1 n in nra 


2.21 


2.18 

n on 1 r» Aoa 


4.90 

n in i n Artd 


1.13 


1.12 

-, -, r , „ ,wfc 



Expt. 0.13 ±0.05" 2.04 ± 0.04" 2.39 ± 0.03 a 3.12±0.02 a 1.04 ± 0.04 b 1.15 ±0.04" 



TABLE VI. AA separation energies, given in units of MeV, of AA-hypernuclei. a) Ref. [E2l 





BAA( A iH) 


Baa(aaH) 


Baa ( A AHe) 


BAA^He) 


Set A, OBE-sim 


0.41 


5.6 


5.5 


14.3 


with CSB 


0.41 


5.2 


6.0 


14.3 


Set A, FSS-sim 


0.53 


6.1 


6.0 


15.1 


with CSB 


0.53 


5.7 


6.5 


15.1 


Set C, OBE-sim 


0.39 


5.5 


5.4 


13.9 



Expt. 10.9 ± 0.6' 



TABLE VII. The spectroscopic factors. Set A AN and OBE-sim AA potentials are used. 





in 


A H (0 + ) 


iH (1 + ) 


iHe 


aaH 


aaH 


AAHe 


0.991 


0.986 


0.992 


0.994 


0.944 


0.962 


0.985 





24 



TABLE VIII. Scattering lengths, in units ol frn, predicted by various AN interactions including CSB. 



Present 
NSC97f|l(J 
NSC89 @| 
NF 01 
ND jS 
Ref. §3 

Ref. i] a ^ N) n > 

O (without U\-(Jn) 



Ap 

-2.83 

-2.51 

-2.73 

-2.18 
-1.77 ±0.28 
-2.16 ±0.26 

-1.83 

-2.45 



An 

-2.26 

-2.68 

-2.86 

-2.40 
-2.03 ±0.32 
-2.67 ±0.35 

-2.45 

-1.83 



Ap 

-1.36 

-1.75 

-1.48 

-1.93 
-2.06 ±0.12 
-1.32 ±0.07 

-1.77 

-1.94 



An 

-1.06 

-1.67 

-1.24 

-1.84 
-1.84 ±0.10 
-1.02 ±0.05 

-1.61 

-1.47 



TABLE IX. A and AA separation energies, given in units of MeV, of A = 5 and 6 hypernuclei as a function of the baryon 
size b (in fm). The case of b = is no quark Pauli calculation. Set A AN and OBE-sim AA potentials are used, a) Ref. [B3J, 
b) Ref. ||. 







Theory 






Expt. 




6 = 


0.6 


0.7 


0.86 




B A ( A He) 


4.98 


4.90 


4.70 


3.64 


3.12±0.02 a 


Baa(a°aH) 
-BaaCaaHc) 

7-> 1 6 TT-\ 


5.6 

5.5 


5.6 
5.5 


5.5 
5.5 


5.1 
5.1 


inn 1 n r-b 



gAA(AAHe) 14.3 13.9 13.0 9.5 10.9 ± O.e 77 



TABLE X. A and AA separation energies (in MeV) for different potentials. The baryon size is set to b = 0.86 fm. 

5a( a Hc) Baa(a^H) Ba_a (aaHb) 

Set A, OBE-sim 3.64 5.1 9.5 

Set C, OBE-sim 3.52 5.0 9.0 

Set A, FSS-sim 3.64 5.5 9.9 



TABLE XL The energy change due to quark Pauli effects. The baryon size is set to b = 0.86 fm. Energy is given in units of 
MeV. 5E denotes the binding energy change B A (b = 0) - B A (b = 0.86) for A He or B AA (b = 0) - Baa (6 = 0.86) for AA H. 







e 2 


<V>f|#|t/>f) 




SE 


iHe 


set A 


0.00406 


513 


2.1 


1.34 


set C 


0.00331 


969 


3.2 


1.37 




set A, OBE-sim 


0.00117 


827 


0.97 


0.5 


AA H 


set A, FSS-sim 


0.00184 


522 


0.96 


0.6 




set C, OBE-sim 


0.000815 


1511 


1.2 


0.5 
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